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Abstract. 

We compare different schemes for the treatment of heavy quark production in Deep- 
Inelastic Scattering (DIS). For fully-integrated quantities such as F 2 {x,Q‘^), we advo¬ 
cate the use of the General-Massive Variable-Flavor-Number (GM-VFN) scheme; we 
present some results showing the progress of a Next-to-Leading Order calculation in 
this scheme. For differential quantities, the Fixed-Flavor-Number (FFN) scheme pro¬ 
vides a more appropriate starting point. We present a new calculation of the azimuthal 
distribution of charm quark production in DIS. All results have been obtained using a 
Monte Carlo program under development. 


Introduction 

The theoretical treatment of heavy quark production has undergone much re¬ 
examination recently. A basic question is whether, or not, the charm quark (or 
bottom quark) should be included in the initial state, with its own parton density 
function (PDF). Traditionally, there have been two approaches used, which we shall 
call the Fixed-Flavor-Number (FFN) scheme and the Zero-Mass Variable-Flavor- 
Number (ZM-VFN) scheme. In the FFN scheme one does not include charm as a 
parton in the nucleon. The charm quark only occurs in the hnal state, and its mass 
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is treated exactly. This is the standard approach for calculations of specihc heavy 
flavor production process, such as charm production in deep inelastic scattering [1] 
and charm or bottom production in hadron-hadron colliders [2], Alternatively, in 
the ZM-VFN scheme one does include a heavy quark PDF in the nucleon at scales 
above the quark mass, but the quark is treated as massless in both the initial and 
hnal states. This is the approach used in most other inclusive jet calculations and 
parton shower Monte Carlo simulations. It also corresponds to the treatment of the 
charm and bottom quark in most global analyses of the parton density functions, 
including DIS structure functions. 

Given these two seemingly contradictory methods of calculating heavy quark 
production cross sections, one may wonder which is correct. The answer, of course, 
is that both are correct at some level. The real question is what is the best method 
for calculating the particular set of observables under consideration at the energy 
scale of the processes being studied. For example, near the charm production 
threshold in deep inelastic scattering, the various energy scales, ~ ~ W^, 

are roughly equal, so that we have essentially a one scale problem. Since the FFN 
scheme treats this scale exactly, including the kinematics of the charm quark mass, 
it should be expected to do well here. However, at higher scales the effect of 
the rric can be neglected, except when it occurs in large logarithms log(Q^/m^). 
These can be resummed by invoking a charm quark parton density function and 
evolving up to Q^. Thus, at higher energies the ZM-VFN scheme will be more 
accurate. As a side remark, we also note that the ZM-VFN scheme is necessary if 
one is to consider the possibility of a nonperturbative charm quark contribution [3] 
to the nucleon PDF beyond that which is generated dynamically by the splitting 
of gluons. 

In general we want a scheme that does the best job at including the largest 
corrections from higher orders in perturbation theory. Since much of the data is 
in an intermediate range of Q^, neither the FFN nor the ZM-VFN schemes are 
completely adequate. For this purpose a General-Massive Variable-Flavor-Number 
(GM-VFN) scheme was devised [4]. (It is sometimes referred to as AGOT, for the 
originators.) In this scheme the charm quark PDF evolves with massless evolution, 
but the heavy quark mass is kept everywhere else in the matrix elements and phase 
space. The leading order process 7 *c ^ c from the ZM-VFN scheme and the leading 
order process 'y*g cc from the FFN scheme are both included in the new scheme. 
However, to keep from double-counting, we must also subtract a term which is the 
convolution of the g ^ q splitting function with the 7 *g ^ q process. This is 
shown diagrammatically in hgure 1. Note that, although the subtraction term may 
appear ad hoc here, it is in fact just the massive equivalent to the Altarelli-Parisi 
subtraction of the 1/e pole at next-to-leading order (NLO) in the ZM-VFN scheme. 

By including the three terms shown in hgure 1 with massless evolution of the 
charm PDF, but full mass-dependence in the coefficient functions, we hnd that the 
calculation in the GM-VFN scheme agrees with the ZM-VFN scheme at large 
and with the FFN scheme near threshold. The three terms in hgure 1 constitute a 
Leading Order (LO) calculation if one considers the charm quark PDF to be 0{as)- 






FIGURE 1. Feynman diagrams for LO GM-VFN scheme. 


This makes sense in the absence of nonpertnrbative charm contribntions since the 
charm qnark PDF arises entirely from splitting of glnons and is conseqnently smaller 
than the glnon PDF. Fnrthermore, this method of as connting leads to rednced 
renormalization/factorization scale (/i) dependence, as we shall see. This scheme 
also allows the inclusion of nonpertnrbative charm contributions to the nucleon 
PDF, which should exist at some level. A recently-proposed prescription by Martin 
et al. (MRRS) [5] is similar in spirit to this. They use the conventional counting of 
as and split the heavy quark mass-dependence into a contribution to the splitting 
functions and a contribution to the coefficient functions. See the talk by Olness at 
this conference for a discussion of this issue [ 6 ]. 

For a more detailed discussion of the theory of heavy quark production, we refer 
to the talk by Tung at this conference [7]. In the remainder of the present talk 
we describe the progress in implementing the GM-VFN scheme at next-to-leading 
order (NLO) in a Monte Carlo program. Currently, we have added the 0{as) 
virtual and real ( 7*0 —>■ gc) contributions to the charm excitation process, with 
the appropriate subtraction. This partial-NLO GM-VFN calculation incorporates 
all of the terms that are in the NLO calculation in the ZM-VFN scheme using the 
conventional a^-counting. The 0{as) virtual and real corrections to the gluon- or 
light quark-initiated processes (and subtractions) that also contribute to the NLO 
GM-VFN scheme have not yet been included. However, we can still gain useful 
insight with the results so far. 


Integrated observables — F2{x, Q‘^) 

As discussed above, the GM-VFN scheme should be the best method of calculat¬ 
ing integrated observables, such as the charm contribution to F 2 {x,Q^). In figures 
2a and 2b we compare the calculation of F 2 {x, Q‘^) in the various schemes. In figure 
2a F 2 {x,Q‘^) is plotted as a function of for x = 0.01. We see that the NLO 
FFN calculation gives a sizeable increase over the LO FFN (tree-level gluon fusion) 
calculation, while the LO GM-VFN scheme is even higher. The partial-NLO GM- 
VFN, which should be most accurate at high Q^, gives a small correction over LO 
GM-VFN. This indicates that the FFN tends to underestimate the cross section 
at high At intermediate we will need the full-NLO GM-VFN calculation, 
which is underway. Note that the GM-VFN scheme is much more economical than 







FIGURE 2. F 2 {x,Q'^) in various schemes, a) versus Q (GeV). b) versus ^/Q. The NLO FFN 
curves are calculated using the code of Harris and Smith [6]. 


the the FFN scheme, since the largest contributions to the cross section are included 
at LO with substantially less work. 

In hgure 2b we show the /i-dependence of F 2 {x,Q^) for x = 0.01 and Q = 10 
GeV. In this plot the NLO FFN calculation only gives a slight improvement in the 
scale dependence over the LO FFN. However, the LO GM-VFN calculation is less 
sensitive to the choice of scale. The partial-NLO GM-VFN calculation actually 
has larger p-dependence than LO GM-VFN, but the remaining terms in the NLO 
GM-VFN calculation should improve this again. In both hgures 2a and 2b the 
wiggles are just due to Monte Garlo statistics. 

Differential observables — (j) distribution 

There are issues that make the GM-VFN scheme somewhat problematic when 
applied to differential distributions. The initial-state charm quark contributions 
include a resummation of collinear radiation diagrams, with the phase space of the 
anticharm quark always integrated out. Thus, the kinematics cannot be treated 
exactly. For instance, the transverse momentum distribution with respect to the 
7 *P-axis is a delta-function, which can only make sense upon integration. In the 
FFN scheme, however, the cc + X hnal state is always treated exactly to a given 
order. Therefore, we expect that it should give a reasonable prediction for the 
differential distributions, as long as we are not too far from threshold. 

In hgure 4 we present a new plot, using our Monte Garlo in the FFN scheme 
at LO, of the distribution of the azimuthal angle 0. This angle is dehned as the 
angle between the lepton (e —e') plane and the hadron (cc) plane in a frame in 
which the 7 * and the proton P are collinear. By symmetry the distribution will be 
of the form ~ A + B cos(j) + C cos 20. In the plot we have made the cuts 

0.01 < y < 0.7, 10 < < 100 GeV, and p\> 2 GeV, where is the transverse 
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FIGURE 3. Distribution in the azimuthal angle (f>. 


momentum with respect to the 7 *P-axis. We have also summed over the charm 
and anticharm contributions, which removes the cos^ term. Although lab frame 
cuts will alter this distribution, they can easily be included in the Monte Carlo 
program. 

We would like to thank Bryan Harris [8] for providing us with the code to produce 
the NLO FFN plots. 
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